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Gravitational waveforms and production could be considerably affected by gravitomagnetic cor- 
rections considered in relativistic theory of orbits. Beside the standard periastron effect of General 
Relativity, new nutation effects come out when corrections are taken into account. Such cor- 
rections emerge as soon as matter-current densities and vector gravitational potentials cannot be 
discarded into dynamics. We study the gravitational waves emitted through the capture, in the 
gravitational field of massive binary systems (e.g. a very massive black hole on which a stellar ob- 
ject is inspiralling) via the quadrupole approximation, considering precession and nutation effects. 
. . . We present a numerical study to obtain the gravitational wave luminosity, the total energy output 

' and the gravitational radiation amplitude. From a crude estimate of the expected number of events 

, towards peculiar targets (e.g. globular clusters) and in particular, the rate of events per year for 

' dense stellar clusters at the Galactic Center (SgrA*), we conclude that this type of capture could 

give signatures to be revealed by interferometric GW antennas, in particular by the forthcoming 
laser interferometer space antenna LISA. 

^ ' PACS numbers: 04.20.Cv,04.50.-|-h,04.80.Nn 
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INTRODUCTION 



Searching for signatures of gravitational waves (GWs) and achieving a suitable classification of emitting sources 



have become two crucial tasks in GW-science. In fact, the today sensitivity levels and theoretical developments 
are leading toward a general picture of GW-phenomena which could not be possible in the previous pioneering era. 
Experimentally, several GW ground-based-laser-interferometer detectors {\Q~^kH z) have been built in the United 
States (LIGO) fH, Europe (VIRGO and GEO) 0, [1| and Japan (TAMA) Q, and are now taking data at designed 
' sensitivities. A laser-interferometer space antenna (LISA) Q (10~^ ~ 10~^Hz) might fly within the next decade. 

From a theoretical point of view, recent years have been characterized by numerous major advances due, essentially, 
to the development of numerical gravity. Concerning the most promising sources to be detected, the GW generation 
problem has improved significantly in relation to the dynamics of binary and multiple systems of compact objects as 
neutron stars and black holes. Besides, the problem of non-geodesic motion of particles in curved spacetime has been 
I developed considering the emission of GWs d, [7| ■ Solving these problems is of considerable importance in order to 
' predict the accurate waveforms of GWs emitted by extreme mass-ratio binaries, which are among the most promising 
' sources for LISA 0. 

!* ■ From a more genuine astrophysical viewpoint, observations towards the central regions of galaxies have detected 
peculiar compact massive objects which are present in almost all observed galaxies. The occurrence of such systems 
has been revealed thanks to the advance in high angular resolution instrumentation for a wide range of electromagnetic 
^ ' wavelengths. Space-telescopes as HST or ground-based telescope, which use adaptive optics, have been extremely 
useful for studying kinematics of galactic internal regions reaching accuracy of mill-pc for the Milky Way and of pc- 
fractions for external galaxies. The main conclusion of all these studies is that the central region of most of galaxies 
is dominated by large compact objects with masses of the order M ~ 10^ IQ'^Mq. In the case of Milky Way, the 
peculiar object in the Sagittarius region (SgrA*) is of the order M ~ 3 x IO^Mq and it is usually addressed as a 
massive black hole (MBH), also if its true physical nature is far to be finally identified [ol. [Tot. 

In any case, a deep link exists between the central MBH and the geometrical, kinematical and dynamical features 
of the host galaxy. In particular, the MBH is correlated with the global shape of galactic spheroid, with the velocity 
dispersion of surrounding stars, with the mean density and the total mass of the host galaxy. Dynamics of stars 
moving around the MBH has a series of interesting characteristics which are of extreme interest for GW detection 
and production. Due to this occurrence, searching for GWs coming from objects interacting with MBHs is a major 
task for GW interferometry from space and ground-based experiments. 

In this paper, we are going to study the evolution of compact binary systems, formed through the capture of a 
moving (stellar) mass m by the gravitational field, whose source is a massive MBH of mass M where m <^ M. One 
expects that small compact objects (1 2OAf0) from the surrounding stellar population will be captured by these 
black holes following many-body scattering interactions at a relatively high rate 0, ^2] ■ It is well known that the 
capture of stellar-mass compact objects by massive MBHs could constitute, potentially, a very important target for 
LISA However, dynamics has to be carefully discussed in order to consider and select all effects coming from 
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standard stellar mass objects inspiralling over MBHs. 

In a previous paper [l5|, we have shown that, in the relativistic weak field approximation, when considering higher 
order corrections to the equations of motion, gravitomagnetic effects in the theory of orbits, can be particularly signif- 
icant, leading also to chaotic behaviors in the transient regime dividing stable from unstable trajectories. Generally, 
such contributions are discarded since they are considered too small. However, in a more accurate analysis, this is 
not true and gravitomagnetic corrections could give peculiar characterization of dynamics. 

In |15|, Newtonian and relativistic theories of orbits have been reviewed considering, in particular, how relativistic 
corrections affect the "classical" orbits [16, 17]. Equations of motion and phase portraits of solutions indicate that, 
beside the standard periastron precession at order c~^, new nutation effects come out at order and it is misleading 
to neglect them. 

According to these effects, orbits remain rather eccentric until the final plunge, and display both extreme relativistic 
perihelion precession and Lense-Thirring [3, HO] precession of the orbital plane due to the spin of MBH, as well 
as orbital decay. In [2l[, it is illustrated how the measured GW- waveforms can effectively map out the spacetime 
geometry close to the MBH. In [13, H^, the classical orbital motion (without relativistic corrections in the motion of 
the binary system) has been studied in the extreme mass ratio limit m ^ M , assuming the stellar system density 
and richness as fundamental parameters. The conclusions have been that i) the GW-waveforms have characterized 
by the orbital motion (in particular, closed or open orbits give rise to very different GW-production and waveform 
shapes); ii) in rich and dense stellar clusters, a large production of GWs can be expected, so that these systems could 
be very interesting for the above mentioned ground-based and space detectors; iii) the amplitudes of the strongest 
GW signals are expected to be roughly an order of magnitude smaller than LISA's instrumental noise. 

In this paper, we investigate the GW emission by binary systems, in the extreme mass ratio limit, by the quadrupole 
approximation, considering orbits affected by both nutation and precession effects, taking into account also gravito- 
magnetic terms in the weak field approximation of the metric. We will see that gravitational waves are emitted with 
a "peculiar" signature related to the orbital features: such a signature may be a "burst" wave- form with a maximum 
in correspondence to the periastron distance or a modulated waveform, according to the orbit stability. Here we face 
this problem discussing in detail the dynamics of such a phenomenon which could greatly improve the statistics of 
possible GW sources. 

Besides, we give estimates of the distributions of these sources and their parameters. It is worth noticing that 
the captures occur when objects, in the dense stellar cusp surrounding a galactic MBH, undergo a close encounter, 
so that the trajectory becomes ti ght enough that orbital decay through emission of GWs dominates the subsequent 
evolution. According to Refs. [11,[2^), for a typical capture, the initial orbital eccentricity is extremely large (typi cally 
1 — e ~ 10~^— 10"'^) and the initial pericenter distance very small (rp 8 — lOOM, where M is the MBH mass j2a |. 
The subsequent orbital evolution may (very roughly) be divided into three stages. In the first and longest stage the 
orbit is extremely eccentric, and GWs are emitted in short "pulses" during pericenter passages. These GW pulses 
slowly remove energy and angular momentum from the system, and the orbit gradually shrinks and circularizes. After 
^ 10'^ — 10^ years (depending on the two masses and the initial eccentricity) the evolution enters its second stage, 
where the orbit is sufficiently circular: the emission can be viewed as continuous. Finally, as the object reaches the 
last stable orbit, the adiabatic inspiral transits to a direct plunge, and the GW signal cuts off. Radiation reaction 
quickly circularizes the orbit over the inspiral phase; however, initial eccentricities are large enough that a substantial 
fraction of captures will maintain high eccentricity until the final plunge. It has been estimated that about half 
of the captures will plunge with eccentricity e > 0.2. While individually-resolvable captures will mostly be detectable 
during the last ~ 1 — 100 yrs of the second stage (depending on the stellar mass m and the MBH mass), radiation 
emitted during the first stage will contribute significantly to the confusion background. As we shall see, the above 
scenario is heavily modified since the gravitomagnetic effects play a crucial role in modify the orbital shapes that are 
far from being simply circular or elliptic and no longer closed. 

The layout of the paper is the following. In Sect. II, we give a summary of gravitomagnetic corrections to the 
metric showing how the geodesic equation results modified by their presence. Besides we study in details orbits with 
such corrections showing the phase portraits and the velocity fields determined by the motion of mass m around 
the MBH. The GW-luminosity in the quadrupole approximation is discussed in Sect. Ill while GW-amplitude with 
gravitomagnetic corrections is discussed in Sect. IV giving also a resume of numerical results. Rate and event number 
estimations are given in Sect.V. Conclusions are drawn in Sect. VI. 

II. GRAVITOMAGNETIC CORRECTIONS 

In a previous paper [lB|, we studied how the relativistic theory of orbits for massive point-like objects is affected by 
gravitomagnetic corrections. In particular, we considered the corrections on the orbits of higher-order terms in v/c 
and this is the main difference with respect to the standard gravitomagnetic effect discussed so far where corrections 
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are taken into account only in the weak field limit and not on the geodesic motion. The problem of gravitomagnetic 
vector potential, entering into the off-diagonal components goi of the metric g^^, can be greatly simplified and the 
corrections can be seen as further powers in the expansion in c^^ (up to c^^). Nevertheless, the effects on the orbit 
behavior are interesting and involve not only the precession at periastron but also nutation corrections. Here we 
briefly recall such previous result sending to for a detailed analysis. 

The metric, in weak field limit where gravitomagnetic corrections are present, is: 



ds^ - + ^) c^dt^ - ^^cdtdx^ - (l - 5) ^ijdx'dx^ , 
where the gravitational Newtonian potential <i>(x) is 



and the vector potential V'' is 



V' = -G f y^^d^x' , (3) 
J |x-x'| 

given by the matter current density pv'' of the moving bodies. This last potential gives rise to the gravitomagnetic 
corrections. 

It is clear that the approximation is up to c^^ in the Taylor expansion. From Eq.([T]), it is straightforward to 
construct a variational principle from which the geodesic equation follows. Then we can derive the orbital equations. 
As standard, we have 



+ r^.i'^i'^ = , (4) 

where dot indicates the differentiation with respect to the affine parameter. In order to put in evidence the gravit- 
omagnetic contributions, let us explicitly calculate the Christoffel symbols at lower orders. By some straightforward 
calculations, one gets 
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9* A. A 



In the approximation which we are going to consider, we are retaining terms up to the orders $/c^ and /c^. 
It is important to point out that we are discarding terms like (<I'/c^)(9$/9a;'^, {V^ / e')d<^ / dx'^ , / c')dV'' / dx^ , 
[V^ / c^)dV^ / dx"^ and of higher orders. Our aim is to show that, in several cases like in tight binary stars, it is not 
correct to discard higher order terms in v / c since physically interesting effects could come out. A vector equation for 
the spatial components of geodesies accounting for the gravitomagnetic effects is ^5|| : 

(Vp 2 4 

= __[V$-e(e.V$)] + -[eA(VAV)] . (6) 

The gravitomagnetic term is the second in Eq.® and it is usually discarded since considered not relevant. This is not 
true if v/c is quite large as in the cases of tight binary systems or point masses approaching to black holes. Orbits, 
corrected by such effects can be explicitly achieved. 
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A. Orbits with gravitomagnetic corrections 

Orbits with gravitomagnetic effects can be obtained starting from the classical theory and then correcting it by 
successive relativistic terms. In [T5| . it is shown that, taking into account the gravitomagnetic terms, in the weak 
field approximation and in the extreme mass-ratio limit m <C M, one obtains a motion with precession and nutation 
by solving numerically the Euler-Lagrange equations. It is possible to obtain the parametric orbital equations of a 
massive particle starting from a variational principle where the canonical Lagrangian is derived by the metric ([1]). 
Being ^ = 0, we have ^ [^] = and then ^ = E where i5 is a constant that can be interpreted as an energy per 
mass unit. Owing to the dependence of the Lagrangian from 9 and 0, we have, in general, 7^ and, furthermore, 
considering the initial conditions — ^ and 6 = we have 9^0. This means, by straightforward calculations, 

■• _ ic^EGMrr - 16 cos G^lVPcj) {ff - IGG^JVP sin cl^irf - 16 cos G'^M^rcjir^ + IGG^M^r sin (^f(/) 
" ^2 + c6r4 - Ac^G^A'P . 

As it is possible to see from Eq.©, planar motion is possible setting the initial condition f = 0, i.e. for the particular 
case of circular orbits; otherwise orbital motions present precession and nutation corrections. Giving explicitly the 
energy first integral 

i = — ( - 4G^(cos6' + sm9(cos(j) + sin0))f + r c^E - iG^i ( {- sin 6* + cos6i(cos0 + sin0))0 (8) 

c-^r — 2cG^ L L V 

+ sin0(cos0 — sin0)(/)) |, 

and the Euler-Lagrange equations (where the energy first integral can be suitably substituted) we get the differential 
system 
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c [rP + G^x) [9^ + sin^ e^^^) p (9) 



cr (rP + 2Gn) 

-AGjii ^(cos6'(coS(^ -f sin0) — sin 6*)^ + sin6'(cos0 — sin0)(/)j r + cG^P — cG^i 



2 (ccot9{rP + 2Gfi) 94>P + r (2G fi esc 9 {sin - cos + cr {rp + Gfi) (pj^ 

P (rc3 + 2Gnc) ' 

c cos 9P {rP + 2Gfi) sin 94>'^ + r (^AG^icos 9{cos 4> + sin (p) - sin 9)i - 2cr [rp + G^i) 9j 

7.2 (^^3 + 2Gfic) 



(10) 



(11) 



Such a system is highly non linear. For its solution, we have to adopt numerical methods. The main characteristics 

of solutions can be evaluated seeing at FiglT] where the stiffness of the system is evident. 

From Eqs.([51)- PT|) . it is clear that additional terms appear with respect to the classical Newtonian motion. Such 

corrections are not independent from and 9. Obviously, these terms become important as soon as velocities approach 

. V . . 

relativistic regimes and the ratio — is quite large. In some physical situations, e.g. in extreme dense globular clusters 

or around the Ga.aoBc Ce„.., Jch a ,at.„ 1 be the ™ge 10- + 10- „t .eshsible at a„ Q and 

references therein). 

In FigUl we have shown the trend of f, (j), 9, as function of t, and 9js[o which is the trend without gravitomagnetic 
correction. It can be seen that this last plot gives deviation from zero that are essentially null, confirming the planarity 
of orbital motions in absence of gravitomagnetic corrections (the differences have at least four orders of magnitude 
between 9Mo{i) and 9{t)). 

To have a further insight of the gravitomagnetic correction relevance on the relativistic orbital motion, we have 

derived a numerical solution with the following parameters and initial conditions: /i = ■ — = 1.4M0, ro = 500^, 

m + M 

E = 0.95, (/>o = 0, 00 = 00 = ^iij'^o and fo = "Tno*-- 1^ Fig-@7 have plotted 9No{t) — ^ without gravitomagnetic 
corrections and 9Grav{t) — f , with gravitomagnetic corrections. In the bottom panel, there is, starting from the left 
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to the right, the trend of the difference between the orbital radii rcrav and rjvo with and without gravitomagnetic 
corrections respectively; we plotted also the differences rQ^av — fNO and tQ^av — iwolred lines) and the ratio between 

coordinated time — versus proper time t (blue line) . It is interesting to see the discrepancy from f of 6* with and 

without the gravitomagnetic effect. It is evident that we have planar orbital motion in the Newtonian case, whilst, in 
presence of gravitomagnetic corrections, there is a tendency to precession and nutation of the orbital plane which give 
rise, orbit by orbit, to cumulative effects (a difference of five orders of magnitude between zj^joit) and zcravit) can be 
evaluated). At the beginning, the effect is very small but, orbit by orbit, it grows and, for a suitable interval of time, 
the effect cannot be neglected [see Fig.([3|), left bottom panel in which the differences in x and y are shown starting 
from the initial orbits up to the last ones by step of about 1500 orbits). On the bottom right, it is shown the basic 
orbit. For about 4850 orbits and a time interval of about 1.7 years, we found that the differences in coordinated time, 
computed with and without gravitomagnetic effects, are increasing as well as the differences in x, y and z coordinates. 
See also Fig. ([¥]) in which wc show the differences between the gravitational wave strain amplitude computed with 
and without the gravitomagnetic orbital corrections (see the discussion in forthcoming Sect. IV). 
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Figure 1: Plots of r = f(t)(left upper panel), (/> = (j){t) (right upper panel), 9 — 9{t){left bottom panel), 6'jvo = dNo{t) (right 
bottom panel). As it is possible to see we have stiff equations owing to the turning points of the orbits. The example we are 
showing has been obtained solving the system for the following parameters and initial conditions: /i — 1.4Mq, ro — 500/i, 



E = 0.95, 



0, 



■rfo and fo 
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III. GRAVITATIONAL WAVE LUMINOSITY IN THE QUADRUPOLE APPROXIMATION 

After the discussion of gravitomagnetic corrections on the orbital motion, let us take into account the problem of 
how GW production and waveforms are affected by such effects. To this purpose, let us start with a short review 
of the quadrupole approximation for the gravitational radiation. This is, in our opinion, the best way to see how 
gravitomagnetic effects correct GW luminosity and waveforms. 

It is well known that the Einstein field equations give a description of how the curvature of space-time is related to 
the energy-momentum distribution. In the weak field approximation, moving massive objects produce gravitational 
waves which propagate in the vacuum with the speed of light. In this approximation, we have 

gfj^^ = + Khf_,^, (|/i^^|<<l), (12) 
being k the gravitational coupling. The field equations are 

Ohf_,^ = -^KTf,^ (13) 



where 
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Figure 2: Plots of 6Mo{t) — § (left upper panel) and of Ocravii) — f (right upper panel). In the bottom panels, rcrav — tmo 
(left) and tcrav — tNO (right) are plotted (red lines). The ratio between coordinated time versus proper time r is also 

plotted (blue line). The examples we are showing have been obtained solving the system for the following parameters and 
initial conditions: fi = IAMq, ro — 500pi, E = 0.95, 0o = 0, 9o — f, 9o — 0, (j>o = —j^ro and ro = ^iJq- 



hi_,i, = hi^,^ - -j^5^uh\x , (14) 

and is the total stress-momentum-energy tensor of the source, including the gravitational stresses. 
A plane GW, solution of (fT^ . can be written as 

h^ii, = h^i, = he^^ cos{u;t - k • x) (15) 

where h is the amplitude, to the frequency, k the wave number and e^i, is a unit polarization tensor, obeying the 
conditions 

Assuming a gauge in which e^i/ is space-like and transverse, a wave travelling in the z direction has two possible 
independent polarizations: 

ei = -^(ii-yy) 62 = -^(xy- yi). (17) 

One can search for wave solutions of p3p . generated by a system of masses undergoing arbitrary motions, and then 
obtain the radiated power. The result, assuming the source dimensions very small with respect to the wavelengths 

dE . . . 

(i.e. the quadrupole approximation J,_7]), is that the power -^^j radiated in a solid angle with polarization e.y, is 

dE _ G fd^Q,, y 

dn " SttcS ^ dt^ 
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Figure 3: Plots of zisio{t) (left upper panel) and zcrav{t) (right upper panel). It is interesting to see the differences of about 
five orders of magnitude between the two plots. At the beginning, the effect is very small but, orbit by orbit, it grows and, 
for a suitable interval of coordinated time, the effect cannot be neglected (see the left bottom panel in which the differences 
in X and y, starting from the initial orbits up to the last ones, by steps of about 1500 orbits, are reported). The internal red 
circle represents the beginning, the middle one is the intermediate situation (green) and the blue one is the final result of the 
correlation between Aa; versus Ay, being Aa; = xorav — xno and Ay — yorav — yNO- On the bottom right, it is shown the 
basic orbit. 
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Figure 4: Plot of the differences of total gravitational waveform h, with and without the gravitomagnetic orbital correction for 
a neutron star of IAMq orbiting around a MBH . The waveform has been computed at the Earth-distance from SgrA* (the 
central Galactic Black Hole). The example we are showing has been obtained solving the systems for the following parameters 
and initial conditions: ^ « IAMq, ro, E = 0.95, ipo = 0, 9o = f , = 0, <;!>o = —jq^o and fo = ~ i^- H is worth noticing that 
frequency modulation gives cumulative effects after suitable long times. 



where Qij is the quadrupole mass tensor 

<3y = - Stjrl) , (19) 

a 
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Figure 5: Plots along the panel lines from left to right of field velocities along the axes of maximum covariances, total 
gravitational emission waveform h and gravitational waveform polarizations h+ and /ix for a neutron star of IAMq. The 
waveform has been computed for the Earth-distance from Sagittarius A (the central Galactic Black Hole SgrA*). The plots 
we are showing have been obtained solving the system for the following parameters and initial conditions: n « IAMq, 
E — 0.95, 4><i = 0, So = f , do = 0, (/!>o = —jq^'o and ro = — i^- From top to bottom of the panels, the orbital radius is 
ro = 20m, ISOO^i, 2500/i. See also Table I. 



ra being the modulus of the vector radius of the a — th particle and the sum running over all masses in the system. 
We must note that the result is independent of the kind of stresses present into dynamics. If one sums Eq. ([T8| over 
the two allowed polarizations, one obtains 

dt^ dt^ ' dt^ dt^ 2 dt^ ) 

where n is the unit vector in the radiation direction. The total radiation rate is obtained by integrating Ea. (|20p over 
all emission directions; the result is 



^ dE 

^dn 



pol 



G 



dE f d^Q,, d^Q,, 1 d^Qu d^Q jj \ 

dn c5 ^ dt^ dt^ 3 dt^ dt^ J ' ^ ' 

It is then possible to estimate the amount of energy emitted in the form of GWs from a system of massive interacting 



Figure 6: Plots along the panel lines from left to right of field velocities along the axes of maximum covariances, total 
gravitational emission waveform h and gravitational waveform polarizations h+ and hx for a Black Hole (BH) of IOMq. The 
waveform has been computed for the Earth-distance to SgrA*. The plots we are showing have been obtained solving the system 
for the following parameters and initial conditions: fj, « IOMq, E — O.95,0o = 0, do ~ f ~ 0,(/f>o = —jq^o and ro = 
From top to bottom of the panels, the orbital radius is ro = 20/i, 1000/^, 2500/^. See also Table I 



objects. In this case, the components of the quadrupole mass tensor are 



Qxx 


= Air2 (3 cos 20 sin 261-1) , 


Qyy 


= /ir^ (3 sin 20 sin 2 6* — 1) , 


Q zz 


= ir2/i(3cos26' + 1) , 


Qxz 


— Qzx ~ r'^ ^J'[^cos4>sv^\29) , 


Qyz 


= Qzy = J'^^d sin26'sin0) , 


Qxy 


= Qyx = r'^^J■ (1 sin^ e sin 20) 



where the masses nii have polar coordinates {risin0cos0, r^sin^sin^, ricos^} and fi is the reduced mass. The 
origin of the motions is taken at the center of mass. Such components can be differentiated with respect to time as in 
Eq. ipTj) . The gravitomagnetic corrections affect, essentially, these quantities and, consequently, the GW amplitude h 
and the radiation rate P as we will see below. 
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IV. GRAVITATIONAL WAVE AMPLITUDE WITH GRAVITOMAGNETIC CORRECTIONS 



Direct signatures of gravitational radiation are given by GW-amplitudes and waveforms. In other words, the 
identification of a GW signal is strictly related to the accurate selection of the waveform shape by interferometers 
or any possible detection tool. Such an achievement could give information on the nature of the GW source, on the 
propagating medium, and, in principle, on the gravitational theory producing such a radiation [27| . 

Considering the formulas of previous Section, the GW-amplitude can be evaluated by 

on 

R being the distance between the source and the observer and, due to the above polarizations, {j, fc} — 1, 2. 

From Ea. (|23p . it is straightforward to show that, for a binary system where m <C M and orbits have gravitomagnetic 
corrections, the Cartesian components of GW-amplitude are 



(3 cos^ cj) SIT? e-l)f'^ + Qr{e cos^ sin 261-0 sin^ 9 sin 20 



+r ( (3 cos^ sin^ 61 - l) f + 3r [6^ cos 261 cos 



smi 



(sin 61 (02 



cos 20 - 



) cos sm ( 



' sin 26* sin 20 



9 cos 9 cos 



(24) 



=2fi (3 sin^ 9 sin^ - l) _^ 6r ^0 sin 20 sin^ 9 + 9 sin 261 sin^ 

+ r ((3 sin^ 9 sin^ - l) f + 3r (j)^ cos 26* sin^ + 0^ sin 261 sin 20+ 
+ sin 9 (j) cos 9 sin^ + sin 9 ^0^ cos 20 + cos sin ( 



(25) 



l^xy ^ ^yx 20 sin 9 (4^0 COS 61 + sin 9^ + 2r ^20 cos 20 sin^ 6I + sin 26* sin 20 

+ i sin 20 (2f sin^ 9 + r{t) (29'^ cos 26^ - 40^ sin^ 9 + 9 sin 26*) ) r + sin^ 9 sin 20 



(26) 

where we are assuming geometrized units. The above formulas have been obtained from Eqs.®, (jlOp . (Ilip The 
gravitomagnetic corrections give rise to signatures on the GW-amplitudes that, in the standard Newtonian orbital 
motion, are not present (see for example [22. .23i]). On the other hand, as discussed in the Introduction, such corrections 
cannot be discarded in peculiar situations as dense stellar clusters or in the vicinity of galaxy central regions. 

Finally, the expected strain amplitude turns out to be /i ~ {h'^x + f^yy + '^^xyV^'^- particular, considering a 
monochromatic GW, we have two independent degrees of freedom which, in the TT gauge, are /i+ = hxx 



xy 



hyx. We are going to evaluate these quantities and results are shown in Figs. [31[H[5l[n] 



hyy and 



A. Numerical results 



Now we have all the ingredients to estimate the effects of gravitomagnetic corrections on the GW-radiation. Cal- 
culations have been performed in geometrized units in order to evaluate better the relative corrections in absence of 
gravitomagnetism. For the numerical simulations, we have assumed the fiducial systems constituted by a m = 1.4M0 
neutron star or m — lOM© massive stellar object orbiting around a MBH M ~ 3 x 10^ M© as SgrA*. In the extreme 

mass-ratio limit, this means that we can consider ^ ~ of about 11 w 1.4Mo and u ~ lOM©. The computations 

m + M 

have been performed starting with orbital radii measured in the mass unit and scaling the distance according to the 
values shown in Table I. As it is possible to see in Table I, starting from rg = 20fj, up to 2500/j,, the orbital eccentricity 

e = — fmm g.^Qjygg towards a circular orbit. In Table I, the GW-frequencies, in mHz, as well as the h amplitude 

strains and the two polarizations /i+ and hx are shown. The values are the mean values of the GW h amplitude 
strains {h = ) and the maxima of the polarization waves (see Figs. [S]and[ni). In Fig. [3 the fiducial LISA 

sensitivity curve is shown [sl] considering the confusion noise produced by White Dwarf binaries (blue curve). We 
show also the h amplitudes (red diamond and green circles for /j, « I.4M0 and ~ IOMq respectively). It is worth 
noticing that, due to very high Signal to Noise Ratio, the binary systems which we are considering result extremely 
interesting, in terms of probability detection, for the LISA interferometer (see Fig. [7]). 
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Table I: GW-amplitudes and frequencies as function of eccentricity e, reduced mass fi, orbital radius ro for the two cases of 
fiducial stellar objects m ~ IAMq and m ~ WMq orbiting around a MBH of mass M ~ 3 X lO''Af0. 
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Figure 7: Plot of estimated mean values of GW-emission in terms of strain h for two binary sources at the Galactic Center SgrA* 
with reduced mass /i « 1.4Mq (red diamonds) and jj, ~ 10MQ(green circles). The blue line is the foreseen LISA sensitivity 
curve. The waveforms have been computed for the Earth-distance to SgrA* . The examples we are showing have been obtained 
solving the systems for the parameters and initial conditions reported in Figs. [5] |6] and in Table I. 



V. EVENT NUMBER ESTIMATIONS TOWARDS SAGITTARIUS A 



At this point, it is important to give some estimates of the number of events where gravitomagnetic effects could be a 
signature for orbital motion and gravitational radiation. From the GW emission point of view, close orbital encounters, 
collisions and tidal interactions have to be dealt on average if we want to investigate the gravitational radiation in a 
dense stellar system. On the other hand, dense stellar regions are the favored target for LISA interferometer [3] so 
it is extremely useful to provide suitable numbers before its launching. 

To this end, it is worth giving the stellar encounter rate producing GWs in astrophysical systems like dense globular 
clusters or the Galactic Center. In general, stars are approximated as point masses. However, in dense regions of 
stellar systems, a star can pass so close to another that they raise tidal forces which dissipate their relative orbital 
kinetic energy and the Newtonian mechanics or the weak field limit of GR cannot be adopted as good approximations. 
In some cases, the loss of energy can be so large that stars form binary (the situation which we have considered here) or 
multiple systems; in other cases, the stars collide and coalesce into a single star; finally stars can exchange gravitational 
interaction in non- returning encounters. 

To investigate and parameterize all these effects, one has to compute the collision time tcoii, where 1/tcoii is the 
collision rate, that is, the average number of physical collisions that a given star suffers per unit time. As a rough 
approximation, one can restrict to stellar clusters in which all stars have the same mass m. 

Let us consider an encounter with initial relative velocity Vp and impact parameter b. The angular momentum per 
unit mass of the reduced particle is i = bvQ. At the distance of closest approach, which we denote by Vcoii^ the radial 
velocity must be zero, and hence the angular momentum is L = TcouVmax, where Vmax is the relative speed at rcoii- 
It is easy to show that ^6|] 
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AGuir, 



coll 



(27) 



If we set rcoii equal to the sum of the radii of the two stars, a coUision wiU occur if the impact parameter is less than 
the value of b, as determined by Eq. ([?7|) . 

The function f{va)d^Va gives the number of stars per unit volume with velocities in the range Va+cPva - The number 
of encounters per unit time with impact parameter less than 6, which are suffered by a given star, is f(ya)(P^a times 
the volume of the annulus with radius b and length vq, that is, 



/(Va)7r6 Vod- Va 



(28) 



where vq = |v — v^l and v is the velocity of the considered star. The quantity in Eq. ipSl) is equal to l/tcoii for a 
star with velocity v: to obtain the mean value of l/tcoii, we average over v by multiplying (|28p by /(v)/j/, where 
V — j f{'v)d^v is the number density of stars and the integration is over d^v. Thus " 



tcoii SttVS 



(29) 



Replacing the variable Va by V = v — Va, the argument of the exponential is then 
we replace the variable v by Vcm = v — -V (the center of mass velocity), then one has 



1 



tcoll STT^crS 

The integral over v^m is given by 



(30) 



(31) 



Thus 



1 7ri/2j. 



tcoll 'icr^ 

The integrals can be easily calculated and then we find 



Tcoii^ +4:GmVr 

coll 

)dV 



(32) 



= A^/nvar^^ii + . (33) 

tcoll cr 

The first term of this result can be derived from the kinetic theory. The rate of interaction is vYi (V), where E is the 
cross-section and (V) is the mean relative speed. Substituting S = i^r'^gn and {V) = Au/^/tt (which is appropriate 
for a Maxwellian distribution whit dispersion a) we recover the first term of (l33l) . The second term represents the 
enhancement in the collision rate by gravitational focusing, that is, the defiection of trajectories by the gravitational 
attraction of the two stars. 

If r* is the stellar radius, we may set rcoii — "^r^. It is convenient to introduce the escape speed from stellar surface, 

= \l and to rewrite Eq. p3|l as 

V 

r = ^ = IQV^varl + - l&^varlil + 9), (34) 



where 



13 



is the Safronov number [iG*]. In evaluating the rate, we are considering only those encounters producing gravitational 
waves, for example, in the LISA range, i.e. between 10"'* and 10^* Hz (see e.g. ^]). Numerically, we have 



3 



If 6 >> 1, the energy dissipated exceeds the relative kinetic energy of the colliding stars, and the stars coalesce into 
a single star. This new star may, in turn, collide and merge with other stars, thereby becoming very massive. As its 
mass increases, the collision time is shorten and then there may be runaway coalescence leading to the formation of 
a few supermassive objects per clusters. If 6 << 1, much of the mass in the colliding stars may be liberated and 
forming new stars or a single supermassive objects (see [ISH^])- Both cases are interesting for LISA purposes. 

Note that when one has the effects of quasi-collisions (where gravitomagnetic effects, in principle, cannot be dis- 
carded) in an encounter of two stars in which the minimal separation is several stellar radii, violent tides will raise on 
the surface of each star. The energy that excites the tides comes from the relative kinetic energy of the stars. This 
effect is important for > > 1 since the loss of small amount of kinetic energy may leave the two stars with negative 
total energy, that is, as a bounded binary system. Successive peri-center passages will dissipates more energy by GW 
radiation, until the binary orbit is nearly circular with a negligible or null GW radiation emission. 

Let us apply these considerations to the Galactic Genter which can be modelled as a system of several compact 
stellar clusters, some of them similar to very compact globular clusters with high emission in X-rays [sij. 

For a typical globular cluster around the Galactic Center, the expected event rate is of the order of 2 x 10"^ yrs~* 
which may be increased at least by a factor ~ 100 if one considers the number of globular clusters in the whole Galaxy. 
If the stellar cluster at the Galactic Genter is taken into account and assuming the total mass M ~ 3 x 10^ M©, the 
velocity dispersion a ~ 150 km s~* and the radius of the object i? ~ 10 pc (where 9 = 4.3), one expects to have 
~ 10"^ open orbit encounters per year. On the other hand, if a cluster with total mass M ~ 10^ M0, a ~ 150 km 
s~* and i? ~ 0.1 pc is considered, an event rate number of the order of unity per year is obtained. These values could 
be realistically achieved by data coming from the forthcoming space interferometer LISA. As a secondary effect, the 
above wave- forms could constitute the " signature" to classify the different stellar encounters thanks to the differences 
of the shapes (see Figs. [S]and|n]). 



VI. CONCLUSIONS 



In this paper we have discussed gravitomagnetic effects on orbital that could give rise to interesting phenomena in 
tight binding systems such as binaries of evolved objects (neutron stars or black holes). The effects become particularly 
relevant when such objects orbit around or fall toward very MBHs as those at the center of galaxies. The effects 
reveal particularly interesting if v/c is in the range (10~* 10~^)c. Gravitomagnetic orbital corrections, after long 
integration time, induce precession and nutation effects capable of affecting the stability basin of the orbits. The 
global structure of such a basin is extremely sensitive to the initial radial velocities and angular velocities, the initial 
energy and masses which can determine possible transitions to chaotic behavior. In principle, GW emission could 
present signatures of gravitomagnetic corrections after suitable integration times in particular for the on-going LISA 
space laser interferometric GW antenna. 
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